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Abstract 

In this work, new types of EWs are introduced. They are constructed by using real 
skew-symmetric operators defined on a single party subsystem of a bipartite d®d system 
and a maximal entangled state in that system. A canonical form for these witnesses is 
proposed which is called canonical EW in corresponding to canonical real skew-symmetric 
operator. Also for each possible partition of the canonical real skew-symmetric operator 
corresponding EW is obtained. The method used for d <8> d case is extended to d\ <8> c?2 
systems. It is shown that there exist distinct possibilities to construct EWs for a 
given di <8> d2 Hilbert space. The optimality and nd-optimality problem is studied for 
each type of EWs. In each step, a large class of quantum PPT states is introduced. It 
is shown that among them there exist entangled PPT states which are detected by the 
constructed witnesses. Also the idea of canonical EWs is extended to obtain other EWs 
with greater PPT entanglement detection power. 

Keywords: Canonical Entanglement Witness, Skew-Symmetric Matrices, En- 
tangled PPT states. 

PACS: 03.67.Mn, 03.65.Ud 
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Entanglement is one of the essential features of quantum physics which has no analogous 
in classical one. Entanglement lies in the heart of the quantum information and quantum 
computation. It is used as a physical resource which allows to realize various quantum infor- 
mation and quantum computation tasks such as quantum cryptography, teleportation, dense 
coding, and key distribution (UE1[3]. The most important problem in quantum entanglement 
is determining boundary of the separable states and entangled ones, which is still not well 
characterized. Although the famous Peres- Horodecki criterion based on positive partial trans- 
pose (PPT) explicitly determines this boundary for low dimensional bipartite systems such 
as 2 Cg) 2 and 2 cg> 3 [U H], but it has no efficiency for entangled PPT states which appear in 
the higher- dimensional compound quantum systems. On the other hand, the entangled PPT 
states belong to the group of entangled states (bound entangled states) which can not be dis- 
tilled. Therefore, distinguishing these states from entangled states that can be distilled (free 
entangled states) is of great importance in quantum communication theory [6]. 

The most general approach to solve and characterize the separability problem in quantum 
theory for any higher-dimensional physical systems is based upon the notion of entanglement 
witness (EW) [7J [8] . The EWs are fundamental tool in entanglement theory since it has been 
shown that, for any entangled state, there exists at least one EW which detects its entanglement 
pO, [9]. A Hermitian operator W is said to be an EW if and only if for all separable states p sep , 
Tr{Wp sep ) > and at least for one entangled state p ent , Tr(Wp ent ) < ( one says that p ent 
is detected by W). Clearly, the construction of EWs is a hard task, although it is easy to 
construct an operator W which has negative expectation value with an entangled state, but it 
is very difficult to check that Tr(Wp) > for all separable states p [TOj [HI [12j [13]. There are 
two types of EWs: decomposable EW (d-EW) that can not detect any entangled PPT states 
and non-decomposable one (nd-EW) which can detect at least one entangled PPT state. It 
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turns out that the development of appropriate separability criteria reduces to construction of 
nd-EWs to detect various entangled PPT states. There have been several considerable efforts 
in constructing nd-EWs (analytically and numerically) for higher- dimensional systems (see e.g. 

[miigngimiiansiEBED). 

In this paper, we describe a method to construct nd-EWs which operate on a bipartite 
Hilbert space with arbitrary dimension. Though Breuer attempted in [18] to build such wit- 
nesses on the basis of time reversal transformation but his work was restricted to those physical 
systems which had angular momentum symmetry. The approach proposed in this paper is gen- 
eral and practical for every physical system. It is based on a skew-symmetric operator which 
acts on a d-dimensional single party subsystem and a maximal entangled state which are called 
characteristic elements of constructed witness. In matrix representation, as will be seen, any 
skew-symmetric matrix U satisfies the relation (a*|U|o:) = (for every state vector \a) that 
lives in the single party subsystem). It is shown that, by using this property of skew-symmetric 
matrices, the related EWs can be constructed analytically for any d <8> d systems. From the 
theory of matrices |22j . the rank of any skew-symmetric matrix is equal to an even number 
(2n). Also every skew-symmetric matrix U can be written as U = QJQ* in which J is the 
canonical form of U (or canonical skew-symmetric matrix) and Q is an orthogonal matrix 
(QQ* = Q*Q = Id)- A witness which is established by J is called canonical EW. It is shown 
that a large number of witnesses are obtained by doing orthogonal transformations specified 
by Q, on the characteristic elements of canonical EW. Thus one is interested to study and 
analyze the canonical EWs in details. 

To develop the idea of the canonical EW for various ranks that J can achieve, the related 
canonical EW is derived. It is proved that for all possible values of rank(J), (rank(J) = 
0, ...In = d), the corresponding canonical EWs is optimal. It is also shown that for full-rank 
J, the related EW is optimal nd-EW and, when the rank(J) < 4, it is optimal d-EW. On 
the other hand, we introduce a new class of PPT states among which there exist entangled 
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PPT states detected by the above mentioned witnesses. By the Jamiolkowski isomorphism [23] 
between operators and maps, the positive maps corresponding to canonical EWs are obtained. 
In the next step, it is assumed that J is full-rank (that is d — 2n). In this case, J is partitioned 
as a direct sum of block-diagonal matrices, which corresponds to a partition of n. As illustrated 
, all of witnesses obtained according to possible partition of n are optimal EWs. Also in this 
step for every partition of n, we construct entangled PPT state detected by the nd-EW which 
is given in the same partition. It is shown that our method can be extended to the bipartite 
d\ ® d2 systems. The shape of canonical EWs are similar to the those ones obtained for d <8> d 
systems. In this case, PPT states are also introduced and by the constructed canonical EWs, 
it is shown that some of them are entangled. 

The paper is organized as follows: In section 2, we investigate our approach to construct 
EWs for d® d systems. The canonical EWs along with the PPT states are introduced. Op- 
timality and nd-optimality for canonical EWs are discussed and also the positive maps corre- 
sponding to canonical EWs are obtained. In section 3, the canonical EWs corresponding to 
possible partition of J and also the related PPT states are discussed. Section 4 is devoted to 
describe the canonical EWs and corresponding PPT states for d\ £g> d2 systems. Also in this 
section, the idea of canonical EW is extended to obtain other EWs. In the end, we summarize 
our result and present our conclusions. 

2 Entanglement Witnesses For d® d Systems 

In this section, we are going to describe a method for constructing EWs. Let us consider H, 
as a d-dimensional Hilbert space devoted to a single particle subsystem. U is defined as a real 
skew-symmetric operator (U T = — U) which acts on the "H. U in matrix form is 



d-i 





(2.1) 
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It is easy to see that for every state \cx) G M. and for every real skew-symmetric matrix U the 
following relation is satisfied 

(a*|U|a) = 0, (2.2) 

i.e., \ct*) and U|a) are orthogonal to each other. Now let us introduce the following Hermitian 
operator 

W = I d ® I d - d\^){^\ - ^(U T ® I,)|^)(^r A (U ® I d ), (2.3) 

where Id <8> Id is the identity operator, T is transposition and \ip), as defined below, is the d® d 
maximal entangled Bell-state 

^ d-l 

l^) = ^Z)l ii )' i = 0,...,d-l. (2.4) 

i=0 

W has the following expectation value with an arbitrary product state \rj) ® \Q) 

( V \ ® (CIW[77> ® |C) = 1 - \((\V*}\ 2 - KC|Uk)| 2 , (2.5) 

in which |t/*) and U|?y) are orthogonal to each other. It is assumed (i]\XJ T XJ\ri) < 1. If UJ77) 
is normalized then \0, in general, can be written as \() = a\rj*) + f3XJ\rj) with \a\ 2 + \(3\ 2 = 1 
otherwise it can not. Consequently, in both cases, the expectation values of the W with respect 
to all separable states are positive, so the Hermitian operator W can be considered as an EW. 
We will show that this type of EWs have ability to detect entangled PPT states, so they are 
nd-EWs. It should be noted that, if we omit the last part of the EW in (12. 3p . the remainder 
part is the reduction EW [TO]; in other words, 

W = W red - d(XJ T ® IdM(^f A (U <g> Id), (2.6) 

where 

W red = Id® Id- 40) (# (2-7) 



Entanglement Witnesses 7 
2.1 Canonical Entanglement Witnesses 

In this subsection, we introduce canonical EW and investigate how other EWs are obtained 
from it. To illustrate, we know that, from the theory of matrices [22], every real skew-symmetric 
matrix on the d-dimensional Hilbert space % can be decomposed as 

U = QJQ T , (2.8) 

where Q is an orthogonal matrix, i.e. QQ T = Q T Q = Id, and J is a block-diagonal one which 
is called the canonical form of U 

J = jo ©jl ©j 2 © ••• ©jn-1, (2.9) 

where 



( o O 

y -A, 0, 



(2.10) 



The scalars Ao, Ax,- ■ ■ ,A n _i, which appear in the J, are invariant factors of U (they are invariant 
under orthogonal transformations). It is clear that the rank of every real skew-symmetric 
matrix U is always an even number (2n); therefore, if the matrix U is full-rank, then 2n = d 
and if d is an odd number, U can not be full-rank anywise. Also it is clear that the eigenvalues 
of U are complete imaginary or zero, but, if U is full-rank, then all of its eigenvalues will be 
imaginary. In addition, the condition (?7|U T U|?7) < 1 yields JJ T = J T J < Id which leads to a 
condition on AjS, i.e. Aj < 1 for % = 0, • • • , n — 1. Now we are ready to introduce canonical 
EW. Consider the following Hermitian operator 

w c = i d ® hi - d^XVI - d(J T © W>(V;f A (J © Id) (2.11) 

By the same prescription proposed for W in equation (I2.3p to be as an EW, Wc can also be an 
EW. It is based on two characteristic elements: J and \ip)(if)\. If we do some transformations 
on these operators, we obtain other EWs. To further illustrate, let be transformed to 
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(Q T <g> Irf) |-0> <-0| (Q <8> Id) then W c becomes as 

= l d ® I d - d(Q T ® I d )|V)(V|(Q ® Id) - d(J T Q T ® I d )|^)(^r A (QJ ® Id) (2.12) 

where Q can be any orthogonal matrix corresponding to any orthogonal transformation in the 
d- dimensional Hilbert space. The subscript \ip) denotes the class of operators obtained by any 
local orthogonal transformation Q <g> Li which is preformed on \ifj){i/j\. Now the calculation of 
the expectation values of over all separable states yields the following equation 

(v\ ® <C|W*to> ® |C> = i - KCIQto*>l 2 - KCIQJM 2 - (2.13) 

By the same argument sketched in equation (12.51) . the above mentioned expectation values are 
always positive. Therefore, the class of operators, W^, can be considered as EWs. The other 
class is obtained by transforming J in to QJQ T as 

Wj = Li <g> l d - d\ip)(ip\ - d(QJ T Q T ® I d )|V)(^r A (QJQ T ® I rf ), (2.14) 

where the subscript J denotes the class of all operators obtained by any orthogonal transfor- 
mation Q which is done on J. By the equation ( 12. 8p . the operator in (12.141) is the EW denoted 
in (12. 3 p . On the other hand, we can see that the class of EWs W^, and W j are related to each 
other by local orthogonal transformation (Q ® Id), i.e. 

W^ = (Q T ®Id)Wj(Q®Id), (2.15) 

therefore they are locally equivalent. These results strongly motivate us to go to study canon- 
ical EW, Wc, in detail. Hence the rest of the paper will be devoted to describe the properties 
of the canonical EW. To this end, we see that the action of J on the basis states of the 
d-dimensional single party Hilbert space H is as 

J|2«) = -\2i + l), J|2« + l) = \2i), z = 0,...,n-l, (2.16) 

and 

J|z) = 0, i = 2n,...,d-l. (2.17) 
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By referring to equation ( 12. lip and using the equations ( 12. 16ft and ( 12. 17ft . Wc is obtained in 
expanded form as 

n-l 

W c = - A?)(|2«, 2% + l)(2i, 2% + 1| + |2* + 1, 2z)(2z + 1, 2z| 

i=0 

-|2*,2i)(2* + l,2* + l| - |2i + 1, 2* + l)<2i, 2*'|) 

n-l 

+ ^ (\2i, 2j) (2t, 2j\ + \2i, 2j + 1) (2i, 2j + 1 1 + \2i + 1, 2j> (2* + 1, 2j| 

+|2i + l,2i + l)(2* + l,2j + l|-|2i,2i)<2j, 2j\ - \2i, 2i)(2j + 1, 2j + 1| 
-|2i + 1, 2* + 1) (2j, 2j| - |2i + 1, 2i + l)(2j + 1, 2j + 1| 
-\A.d2i + l,2j)(2j + 1, 2i| - \2i + 1, 2j + l)(2j, 2*| 
-|2i, 2j)(2j + 1, 2i + 1| + |2i, 2j + l)(2j, 2* + 1|)) 

n-l d-1 

+ SZ)(l 2i '^< 2i '^'l + |2i + l,i)<2i + l,j| + |i,2i)(j,2*| + |.7,2* + l)(j,2z + 1| 
-|2*,2*)<j, j| - |2* + 1,2* + - \j,j)(2i,2i\ - \j,j)(2i + 1,2* + 1|) 

+ E(l^'>^'|-|M>0-,J|). (2.18) 

i,j=2n 

Consider the following operator 

n-l 

Ta = ^(1 - Aj)(|2*,2* + 1)(2*,2*+ 1| + |2* + l,2z)(2z + l,2z'| 

i=0 

- |2*,2*)(2* + 1,2* + 1| - |2*+l,2z + l)(2z,2z|). (2.19) 

The Wc is composed of the 0* A and the remainder one which is called Wc(A , Ai, A n _i), 
i.e. 

W c = Ta + Wc(A , Ai, A n _0. (2.20) 

We see that the operator O is positive; therefore, from [23] for every PPT state p, Tr(0*- 4 p) > 0. 
It becomes clear that the following inequality is satisfied for every PPT state 

Tr(Wcp) > Tr(Wc(A , A 1; X n -x)p). (2.21) 
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This inequality enables us to take A; = 1 for % — 0, n — 1 so Ta becomes zero and Wc = 
Wc(l, 1, 1). From now on, we keep discussing on Wc for Aj = 1 (i — 0,...,n — 1) and 
rewriting Wc as 

n-l 

W C =J2 d 2 *' 2 -?>^' 2 ^'l + I 2 *' ^ + X > < 2 ^ 2 1 + 1 \ + I 2 * + !> 2 ^') ( 2i + !> ^1 

+|2i + l,2j + l)(2i + l,2j + l|-|2i,2i)(2^2j|-|2i,2i)(2j + l,2j + l| 
-|2i + 1, 2i + l)(2j, 2j| - |2i + 1, 22 + l)(2j + 1, 2j + 1| 
-|2i + 1, 2j) (2j + 1, 2*| + |2i + 1, 2j + l)(2j, 2i| 
+ |2z,2j>(2j + l,2i + 1| - |2i,2j + l><2j,2i + 1|) 

n-l d-1 

+ EE(l 2i '^< 2i ^'l + |2i + 1, j>(2i + + |j,2i)0 - ,2i| + |j,2z + 1)0', 2i + 1| 
-|2z,2i)(j,j| - |2i + l,2z + l)(j,j| - |j, j><2i,2i| - \j,j)(2i + l,2i + l\) 

+ (2.22) 

i,j=2n 

To disambiguate, W c can be briefly written as 

W c = Of- 4 + T 2 A + Wopc, (2.23) 

where 

n-l d-1 

°? = EEd 2 ^'H 2 ^1 + I 2 * + U>(2i + M + b', 2z><j, 2z| + | j, 2i + 1>0", 2i + 1| 

i=0 j = 2n 

- |2i,2i)(j,j| - |2i + l,2i + l)(j,j| - \j,j)(2i,2i\ - \j,j)(2i + l,2i + l\), (2.24) 

<i-l 

= £(M<M'|-|M>0',.7'I) (2.25) 

i,j=2n 

and 

n-l 

Wopc = ^ (|2i, 2j)(2i, 2j| + |2i, 2j + 1) (2i, 2j + l\ + \2i + 1,2 j) (2i + l, 2j\ 
+\2i + l,2j + l)(2i + l,2j + l\-\2i,2i)(2j,2j\-\2i,2i)(2j + l,2j + l\ 
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-\2i + l,2i + l)(2j, 2j\ - \2i + 1, 2i + l)(2j + 1, 2j + 1| 
-|2i + 1, 2j> (2j + 1, 2i\ + \2i + I, 2j + 1) (2j, 2i| 

+ |2z, 2j) (2j + 1, 2i + 1 1 - \2i, 2j + 1) (2j, 2i + l\). (2.26) 

Note that if the rank of J (2n) is two then Wopc will be zero and therefore by considering 
the equation (12.23 p . it is concluded that Wc is a d-EW. Since we are interested in dealing 
with nd-EWs then those d® d quantum systems for which rank(J) > 4 are discussed. Hence 
we conclude that, by this approach, nd-EWs can be constructed only for those systems with 
d > 4. Clearly if J is full-rank (2n = d) then in equation (12.231) . 0\ A and 0^ A will be zero. 
Now we claim that the Wc* type witnesses are able to detect entangled PPT states so they are 
nd-EWs. 

2.2 PPT states 

This subsection is devoted to construct PPT states and determine a subset of them as a set of 
entangled PPT states whose entanglement are detected by the EWs introduced in the previous 
subsection. Let us write the following operator 

^ n— 1 

p= —(ool^X^I + a ^(|2z,2i)(2i,2z| + \2i + 1, 2% + l)(2i + 1, 2% + 1| 

i=0 

-\2i,2i)(2i + l,2i + 1| - \2i + l,2i + l)(2z,2z|) 

n-1 

+ 5^(a 2i+2! 2i|2i + 2,2z)(2z + 2,2z| +a 2mj2i+3 |2i + l,2i + 3)(2i + l,2z + 3| 

i=0 

-Ci(\2i + 2, 2i) (2i + 1, 2% + 3| + \2i + 1, 2i + 3) (2i + 2, 2i\)) 

n— 1 n— 1 

+ a 2h2j \2i,2j)(2i,2j\+ a 2i+1>2j+1 \2i+l,2j + l)(2i + l,2j + l\ 

i^j=0,i-jj£l,j-i^n-l i^j=0,j-i^l,i-j^n-l 
n-1 

+ ^K, 2i+1 |2i,2j + l)(2z,2j + 1| + a 2l+lj2i |2z + l,2j)(2i + l,2j\) 

i,j=0 

d-l n-1 

+ Yl 5Z(°2^l 2 ^ «>( 2 i> i|+a2i+M|2j+l, i)(2j+l, i|+a ij2j |z, 2j)(z, 2j|+a ii2i+1 |z, 2j+l)(z, 2j+l|) 

i=2n j=0 
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d-l 

+ <kj\i,3){i,3\), (2-27) 

where Af is the normalization factor and equals to 

n— 1 n— 1 

AT = (d + 2n)a + ^ (a 2l _ 2j + a 2i+12j+1 ) + ( 

d-l n-1 d-l 

)+ °W ( 2 - 28 ) 

i=2n j=0 i=ij=2n 

The positivity conditions impose that all of the multipliers which appear in the p are positive 
semi definite and in addition the following inequality must be satisfied 

« 2 * + i, 2!+ 3 a 2 I+ 2, 2l > c l i = 0, ra - 1, (2.29) 

where addition in the subscripts is done by module (2n). Also the PPT conditions are as 
follows: 

a 2 i,2j a 2j ,2i > al, i,j = 0, ...,n - 1, i ^ j, 
OM+i,2j+i a 2 j+i,2i+i > al, i,j = 0, ...,n- 1, z 7^ j, 
a2i,2j+i a 2j +i,2i > al, i,j = 0, ...,n - 1, i ^ j, 

Cl2i+l,2i 0<2i+2,2i+3 > Cf , Z = 0, 77. — 1, 

a 2 j,i a i: 2j > a 2 Q , j = 0, ...,n-l, z = 2n, d - 1, 

«2j+i,i a it2 j+i > al, j = 0,...,n-l, z = 2n, d - 1, 

Oij a j:i > al, i,j = 2n, ...,d- 1, i ^ j. (2.30) 

Therefore, by these two groups of inequalities, the operator p becomes as a density operator 
with positive partial transpose. In the next step, it is shown that the expectation value of 
the witness W c with respect to the p, under positivity and PPT conditions, really fulfills the 
following inequality (the proving of the following Lower bound is given in the appendix A.) 

— 2n 

Tr(VW) > - 5 ; rn r, 71 = 2, 3, 4, • • • . (2.31) 

v Cl ! ~ d + 4n 2 + (d-2n)(d + 2n- 1) ' ' ' ' v ; 
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Hence if the PPT state p satisfies the following inequality 

-2n 



<Tr(W c p)<0, n = 2,3,4,---, (2.32) 



d + An 2 + (d - 2n)(d + 2n-l) 
then it will be an entangled PPT state whose entanglement is detected by Wo- On the other 
hand, the witness Wc which can detect the entangled PPT state p, becomes as a nd-EW which 
proves our claim. It is also seen that the lower bound of the inequality in (I2.27P depends on 
the rank of the matrix J and the dimension of the Hilbert space of the single party subsystem. 
It is also obvious that if J is full-rank, then the lower bound becomes smaller. Finally; when 
the state p violate the PPT conditions, the expectation value of the entanglement witness Wc 
with respect to the density operator p satisfies the following inequality 

4n (n - 1) + (d - 2n)(d + 2n - 1) < ^ , < "2n 



d + 2n ~ v rj d + An 2 + {d - 2n)(d + 2n - 1)' 

ra = 2,3,4,- •- . (2.33) 

.Clearly the lower bound becomes greater when the matrix J is full-rank. 
2.3 Optimal Canonical Entanglement Witnesses 

Now we discuss at first the optimality of canonical EWs by investigating optimality of canoni- 
cal EW when J is full-rank (2n=d). The proving of optimality for the others (2n < d) is similar 
to this one. Secondly, we describe nd-optimality of canonical EWs. There exist different defi- 
nitions of optimal entanglement witness. Our description is based on the definition introduced 
by Lewenstein et.al., [21]. One has two EWs Wi and W2, W2 is finer than Wi if they differ by 
a positive operator P. We say that W is optimal iff for all P and e > 0, W' = (1 + e)W — eP is 
not an EW. P is positive operator and PP\y — where P\y = {It), Tr(W|7)(7| = 0)} in which 
1 7) is separable state. Since any positive operator can be written as a convex combination of 
pure product states so let us assume that P = \if))(ip\ in which 

d-l 
i,j=0 
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By referring to the equation (12. 5 j) and using J instead of U , a typical separable state (7) G P\y c 
has the following form 

\l) = \v)®(a\vl+PJ\v)), (2-35) 

where |cf + |/3| 2 = 1. We define |r?> = ££ffaj2i> + + 1)) such that EIToCKI 2 + 

l^+il 2 ) = 1- For simplicity we choose a = 1 and /3 = (I7) = \rj) <8> |r/*)). Therefore, the 
separable state (7) in expanded form is written as 

n-l 

It) = E ^ 2 ^'> + i»C. I 2 '' ^ + x > 

i,i=0 

+ r) 2i+lV * 2 . \2i + 1, 2j> + |2i + 1, 2j + 1). (2.36) 

It is proven P=0 as: 

1. r} 2 . — 5ik which gives the following separable state 

| 7 > = \2k, 2k) — »- ( 7 |<0> = a 2K2k = 0. (2.37) 

| 7 > = \2k + 1, 2k + 1} — ► (7|V) = a 2fe+1)2fe+1 = 0. (2.38) 
3- 1^ = 0', |a'| 2 + OT = l, 

| 7 ) = \a'\ 2 \2k, 2k) + a'/3'*|2A;, 2k + 1) + /9V|2Jfe + 1, 2Jfe) + |/3'| 2 |2A; + 1, 2k + 1), 

( 7 |^) = /3'a'*a 2 fe + i,2/fc + a'/3'*a2fc,2fc+i = 

— >■ ft2fc,2fc+i = 0, a2fe+i,2fe = 0. (2.39) 
4. r? 2fc = a', r? 2i =/?', k<l, \af + \(3f = l, 

1 7) = |a'| 2 |2A;, 2Jfc) + a'£'*|2A;, 2/) + /3V*|2/, 2fc) + |/3'| 2 |2Z, 2/), 
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( 7 |^) = a'/3'*a 2 k,2i + P' a* 0,21,2k = 

— ► a 2k ,2i = 0, a 2l , 2k = 0. (2.40) 

5. n 2k+1 = a', ?7 2i+1 = p', k<l , |a'| 2 + |/?'| 2 = 1, 

|7) = \a'\ 2 \2k + 1, 2A; + 1) + |2fc + 1, 21 + 1) + /3V|2Z + 1, 2k + 1) + |/3'| 2 |2Z + 1, 2/ + 1), 

(7|V>) = a / /3'*a 2 fc + i,2«+i + /3 / at / *a 2 j+i ) 2fc+i = 

— ► G2fc+i,2i+i = 0, a>2i+i,2k+i = 0- (2-41) 

6. n 2k = a', n 2l+1 = p' , k<l, \a'\ 2 + |/3'| 2 = 1 

| 7 ) = \a'\ 2 \2k, 2k) + a'P'*\2k, 21 + 1) + 0'a/*\2l + 1, 2Jfc) + |/3'| 2 |2/ + 1, 2/ + 1), 

(7|V») = a'p'*a 2 k,2i+i + P'ot'*a 2 i+i ) 2k = 

— ► a 2 fc,2i+i = 0, a 2 i+i,2k = 0. (2.42) 

7. n 2k+1 = a', ^ = p' , k < I, \af + |/3'| 2 = 1, 

| 7 ) = \Pf\2l, 21) + f3'a'*\2l, 2k + 1) + a'p'*\2k + 1, 21) + |a'| 2 |2A; + 1, 2k + 1), 

( 7 |*0) = a'p'*a 2 k+i,2i + (3'a'*a 2 i t2k+1 = 

— > «2/,2fc+i = 0, a 2 k+i,2i = 0. (2.43) 

These equations explicitly show that P=0 and therefore the canonical EW is optimal. The 
proving of optimality for Wc, when J is not full-rank (and specially is zero), is similar to 
the previous one except by noting that the typical separable state (7) G Pc has the form 



|7> = \v) ® \V*) where \n) = Y^=d(vJ 2i ) + V 2l+ i\ 2i + !)) + E;=2n^K)- Therefore, it i 



is 



concluded that the canonical EW (12. lip is optimal for all ranks of J (rank(J) = 0, 2n = d), 
in other words, its optimality is independent from the rank of J. On the other hand, to discuss 
nd-optimality, we define d^ c = {p > 0\p TA > 0, Tr(W c p) < 0}, i.e., the set of entangled 
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PPT states detected by Wc- Given two nd-EWs Wei and Wc2, we say that Wc2 is nd-finer 
than Wei, if ^w C i — ^w C 2) i- e -> if a ^ °f the entangled PPT states detected by Wei are also 
detected by Wc2- We say that Wc is optimal nd-EW, if there exists no other nd-EW which is 
nd-finer than it. So by keeping this in mind, we determine optimal nd-EW among canonical 
EWs. By referring to equation (|2.23p . we see that the operators 1 and 2 are positive, so 
0^ A and Oj A are optimal d-EWs j2l] and Wopc is a canonical EW on the 2n <g) 2n Hilbert 
space which is a subspace of d®d one (the proving that Wopc is an EW, is given in appendix 
B). We know that for any PPT state the following inequality is satisfied 



From the Lewenstein definition of an optimal nd-EW, it is clear that each entangled PPT state 
detected by Wc, is also detected by Wopc- It is proven that Wopc is an optimal nd-EW. 
To this aim we say that a nd-EW, W, is optimal nd-EW iff for all decomposable operator D 
(D = P + Q* A where P and Q are positive operators) and e > 0, W = (1 + e)W — eD is not an 
EW. In the proving of optimality for canonical EWs, It was shown that P = 0. To illustrate 
that Wopc is an optimal nd-EW, it must be shown that for Q Ta P Wopc = then Q Ta = 0. Let 
us assume that Q = \(p)(<p\ in which \ip) = X^7=o a »j'K) j)- As previously, since J is full-rank on 
2n-dimensional subspace of the d-dimensional one party Hilbert space H, a typical separable 
state 1 7) G P\v pc which lies in the 2n £g> 2n subspace (see appendix B), is written as |#) = 
\rf) ® (a\n*) + pj\r))). On the other hand, since Tr(Q* A |a) (a\ <8> |/3) = Tr(Q>*) (a*\ ® |/3) (/3\) 
then we calculate the expectation values of Q with the products = \rj*) <g> (ot\rf) + /3J\rj)). 
I'd) in expanded form is 



Tr(Wcp) >Tr(W Pcp)- 



(2.44) 



n-l 



i,j=0 



(2.45) 
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One can show Q = on the subspace 2n <g> 2n as below 

1. r) 2 . — 5 ik , ol — 1, /3 — which gives the following separable state 

|0) = |2fc, 2A;) — >• (&\<p) = a 2k , 2k = 0. (2.46) 

2. r? 2l = , a = 0, /? = 1 

\#) = \2k, 2k + 1) — ► (%) = a 2fc , 2fc+1 = 0. (2.47) 

3. r? 2l+1 = , a = 0, £ = 1 

= |2fc + 1, 2fc) — > (%) = a 2fc+1 , 2fc = 0. (2.48) 

4. r} 2i+1 = 5 ik , a = = 

|t?> = |2fc + 1, 2fc + 1) — >■ (%) = a 2fe+1)2fe+1 = 0. (2.49) 

5. %k = a',n 2l = f3> , \af + \(3f = 1 , k < I , a = 0,f3 = 1 

| t?) = |a'| 2 |2A;, 2A; + 1) + a'f3'*\2l, 2k + 1) + /5W|2A;, 2/ + 1) + |/3'| 2 |2/, 2/ + 1), 

= a'(3'*a 2lt2k+1 + f3'a'*a 2kt2l+1 = 
— > a 2 i, 2k +i = , a 2fci2 i+i = 0. (2.50) 

6. %k+1 = a',n 2l+1 = P , \af + \f3f = 1 , k < I , a = 0, f3 = 1 

1 7?) = |a'| 2 |2A; + 1, 2fc) + |2/ + 1, 2ifc) + /3W|2A; + 1, 2/) + \f3'\ 2 \2l + 1, 2/), 

(0\<p) = a'P'*a 2l+lt2k + /?'«'* a 2k+h2l = 

— >■ a 2 i+i,2k = 0, a 2k+ i,2i = 0. (2-51) 

7. 77 2fc = a', n 2l+1 = (3> , \af + |/?'| 2 = 1 , k < I , a = 0, (3 = 1 

|#) = -|a'| 2 |2A;, 2Jfe + 1) - a'/?'*|2/ + 1, 2A; + 1) + /3'a / *|2fc, 2/) + |/?'| 2 |2/ + 1, 21), 
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(1%) = -a 'P'*a 2 i + i,2k + /3'a'*a 2k ,2i = 

— > 02i+i,2Jfe+i = 0. a 2 fc,2« = 0. (2.52) 
8. n 2k+1 = a', ^ = 0' , |«'| 2 + \(3f = 1 , k < I , a = 0,0 = 1 

|#) = \a'\ 2 \2k + 1, 2k) + a'P'*\2l, 2k) - 0'a'*\2k + 1, 21 + 1) - |/?'| 2 |2/, 21 + 1), 
(?%) = a'/3'*a 2 i,2k - /3 / a / *a2fe+i,2«+i = 

> «2fc+l,2/+l = 0; &2l,2k = 0. (2.53) 

therefore, Q = on the 2n®2n subspace. From the other side, since all of the separable states 
which lie on the complement subspace of the 2n g) 2n one belong to the P\v pc ( see appendix 
B), then it is obvious that Q is also zero on that subspace. Consequently Q = on d <8> d 
Hilbert space hence Wopc is optimal nd-EW. Clearly, if J is full-rank, then Op 1 and 0^ A will 
be zero so Wc = Wopc and Wopc is an optimal nd-EW on the d <S> d Hilbert space. It is 
concluded that when J is full-rank, Wc is optimal nd-EW and when J is not full-rank Wc 
is not optimal nd-EW; therefore, despite the optimality, the nd-optimality of canonical EWs 
depends on the rank of J. 

2.4 The Positive Maps Corresponding to Canonical EWs 

Since the nd-EWs have an essential role in the studying of separability problem in quan- 
tum theory, by using Jamiolkowski isomorphism [23] between operators and maps, the non- 
decomposable positive maps (or nd-positive maps) have the same role as nd-EWs. By this iso- 
morphism, one can obtain the corresponding positive map of the canonical EW Wc € Hd^Hd 
(12. lip as discussed in subsection (2.1). Consider the following equation 

0(p) = Tr B (Wc(I d ®p T )). (2.54) 

Where p is a density operator on the d-dimensional Hilbert space. This equation shows how to 
construct the map from a given operator Wc- After some calculations, the following result 
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is obtained 



0(p) = I d Tr(p)-p-jVJ- 



(2.55) 



From the properties of Wc discussed earlier, we expect that, if the rank(J) < 4, then the <fi(p) 
is decomposable positive map (or d-positive map), especially when rank(J) = 0, 0(p) is the 
well-known reduction map |25J. Therefore, for rank(J) > 4, <p(p) is nd-positive map. 

3 Entanglement Witnesses Corresponding to possible 
partitions of J 

In this section, by referring to each possible partition of J, we are going to construct a new 
set of canonical EWs. Therefore, for a given J, we have a set of canonical EWs corresponding 
to the set of possible partitions of J. It is shown that, for a given partition of J, a PPT 
state is constructed for that partition. The entanglement of this PPT state in some range of 
parameters is detected by the corresponding canonical EW established in the same partition. 
Suppose that J is full-rank, i.e., d = In. Consider a partition of d-dimensional single party 
Hilbert space % to its 2 /^-dimensional subspaces, %2 W s, through the following direct sum 



Wo = Id ® Id - - <UT ® IdM(V>F A (Ui © I d ) - d{\J T 2 <g> I d M(^F A (U 2 © h) 

- d(X]T ® I d )^)<Vf A (U 3 ® I d ) d(XJl ® IM)(iJj\ t *(XJ v ® I,), (3.57) 



in which every Uj is block-diagonal full-rank matrix on (k = 1, • • • , v) such that 



H — n-211,1 © %/Lt 2 © ^2/i 3 © • • • © 'H.211U- 



(3.56) 



Also consider the following Hermitian operator 



J = Ui + u 2 + u 3 + • • ■ + u 



(3.58) 
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where UjU f = (i ^ j = 1, v). Each Uj is a full-rank canonical skew-symmetric matrix in 
"H 2Ati . Therefore the rank of J is the sum of the ranks of UjS that is 

2n = 2/i! + 2fi 2 + 2// 3 + • • • + 2//„ (3.59) 

hence we obtain the next result 

n = Hi + fi 2 + ^3 H H /V (3.60) 

It is well-known that the numbers (/xi, /x 2 , 1^3, • ■ • , A*i/) are a partition of n. Generally from [26], 
for a given number n there are p(n) number of partitions (/xi, /x 2 , fJ-3, • • • ,(J> U ) with /ix > /x 2 > 
^3 • ' ' > > (f° r example consider 5, the number of partitions for it, is p(5) = 7). Hence 
we say that the U&S in ( 13.581) are a partition of J. Therefore for the other possible partitions 
of n, we have corresponding partitions for J and corresponding Hermitian operators such as 
Wc which is renamed as Wc(/xi, /x 2 , /X3, • • • , (h>)- 

In the same way, as described in section 2, the expectation values of Hermitian operator 
Wc with all product states are given as 

(V\ ® (C|W (Mi,A*3, " " " ,^-i,^)\v) ® 10 = 1 - l(Cl^)| 2 - KClUilr/)! 2 - |(C|U 2 |r/)| 2 

-KClUal^)! 2 |<C|U,|?7>| 2 , (3.61) 

in which the states {\r}*) , Ui|t7) , U 2 |r?) , U 3 1 77) , ••• , XJ u \rj)} are orthogonal to each other. 
If Uj|77), for each i, is normalized then \Q can be written as \Q = o\rf) + Yli=o A^il 7 ?) with 
I or 1 2 + Ym=i IAI 2 = 1 otherwise it can not. Therefore Wc(/xi, /x 2 , /X3, • ■ ■ , /x y ) has positive 
expectation values with all separable states so it can be considered as a canonical EW for the 
partition (/xi, /x 2 , /x 3 , • • ■ , (j, v ). par The action of Uj (i = on the basis states of H 2ai . 

(z = 1, i/) is as 

Ui|2fc) = -|2Jfe + l), Ui|2A; + l) = \2k), 

i = !,...,!/, /c = 0, fii — 1, (3.62) 
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therefore Wc(fa, fa, fa, • ■ ■ , fa) in equation (13. 5 7ft for a given partition (/ii, fa, fa,- ■ ■ , fa) can 
be obtained as 

WcO"i, A*2, A*3, • • • , A*i/) = WcOui, /u 2 , a*3, • • • , fa) 

+ O tA {fa) + 0^(/i 2 ) + 0^(// 3 ) + • • • + O tA (ji v -x), (3.63) 

where 

W c (/ii, fa, fa,--- , fa) = W c (mi) © Wc0 2 ) © Wc(fa) © • ■ ■ © Wc(fa) (3.64) 

The operators Wc(/Uj) (i = 1, • • • , v) and 0* A (^) (z = 1, • • • , v — 1) have been given in the 
appendix C. By the same prescription shown in section (2.3), the Wc(fa, fa, fa, • ■ ■ , fa) is 
optimal EW for each partition of n. Each Wc(/ii) (i — 1, ■ • • , v) is an EW in the 2/Xj ® 2/Zj 
subspace of the d <S> d Hilbert space. Also for every fi { > 1 the corresponding Wc(/ii) is an 
optimal nd-EW and hence W c (fa, fa, fa, - 4 4 , fa) which is the direct some of the Wc(/U«)s, is 
also an optimal nd-EW. The operators 0(/ij) (i = 1, u — 1) are positive operators then the 
operators 0* A (/Zj) (i = 1, ■ ■ • , v — 1) are optimal d-EWs. 

As mentioned in section (2.1), if a given ^ be equal to one then the corresponding Wc(a*0 
will be zero. Therefore if all of /ijS become one, i.e. for the partition (1, 1, 1, ■ • • ,1), then 
Wc(l, 1, 1, ■ • • ,1) will be an optimal d-EW. The witness corresponding to the partition (n = 
fa), i.e. Wc{fa), is an optimal nd-EW which was discussed earlier. In the end, the witnesses 
which correspond to the other partitions between these two partitions , as the equation (13.631) , 
are a mixture of optimal d-EWs and the optimal nd-EWs. Therefore, by considering the 
optimality of the nd-EWs discussed in subsection (2.3), these witnesses are not optimal nd- 
EWs. 

In the next step, for a given partition (fa, fa, fa, ■ ■ - fa), a PPT state is introduced. This 
state can be entangled and the signature of entanglement for it is shown by the witness which 
corresponds to the same partition discussed above. This state is the following one 

p{fa,fa,fa, ■■■ ,fa) 
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(q(Hi,H2, A*3, • • • , to) + cr(A*i) + cr(A*i) + cr(/x 3 ) H h ct^-i)) (3.65) 



N{to,to,to, ■■■ ,to) 
such that 



/i 2 , A*3, • • • , to) = pM © pM © p(^) © • • • © pM, (3.66) 

where M{to, /i2, A*3> • • • , ^) is the norm of the p(p 1 , p 2 , to, ' ' ' ,to)- The operators p(to) s 
(i = ,v) are also (unnormalized) PPT states in the subspace 2^j ® 2//j and by the 

witnesses Wc(/ii)s, they are entangled PPT states. Also the operator q(h±, to, to, ' ' ' 
which is the direct sum of the p(pi)s, is PPT state so its entanglement is detected by the 
witness W c (pi, to, to, ' ' • , to) an d finally cr(//j)s (z = 1, • • • ,v — 1) are Hermitian operators. 
All these operators together with the positivity and PPT conditions for p(to, to, to, ' ' ' , to) 
have been given in appendix D. 

Now at the end of this section the expectation value of Wc(to, to, to, " " ' , to) with respect 
to the PPT state p(pi, to, to, ' " " ,to), by considering both positivity and PPT conditions, 
is calculated. In the same way as in subsection (2.2), after some calculations we obtain the 
following lower bound 

Tr[W c (^i, fi2, to, ■ ■ ■ , to)p(to,to, to,'-' , to)} 

- ~ EH^s; + i) + Ap e {d - 2 ELi io5 ' (3 ' 67) 

Therefore we say that if a PPT state p(pi, to, to, ' ' ' , to) satisfy the following inequality 

d 

"SSS + 1) + 4to(d ~ 2 ELi to)) 

< Tr[W c (to, to, to,-- - , to)p(to,to, to,-- - , to)] < 0, (3.68) 
then it is an entangled PPT state detected by the witness Wc(to, to, to, ' ' ' ,to)- Finally if 
p(to, to, to, ' ' ' , to) violate the PPT condition, it always satisfies the next inequality 
T,:=i(2to(to - 1) + 2to(d - 2 J2l=i to)) 



d 



< Tr[W c (to, to, to,-- - , to)p(to, to, to,-- - , to)} 

< " TUVtoWe + 1) + 4^(d - 2 ELi to)) ' (3 ' 69) 
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In this section, at first, we extend our approach used for d® d cases in the previous sections to 
obtain EWs in d\ <g> d 2 quantum systems [d\ < ci 2 ). Secondly, we extend the idea of canonical 
EW introduced in subsection (2.1) to construct other EWs. Let us introduce a projection 
operator P c on the (i 2 -dimensional party in which Im(P c ) is a dx-dimensional Hilbert space. 
Since in the <i 2 - dimensional vector space one can construct di-dimensional subspace by Cjj 
(the number of combinations of <i 2 distinct objects taken d\ at a time without repetitions) 
distinct ways; therefore, the subscript c (c = 1, • • • , C^) dentes a projection operator among 
the set of Cj ones . From the linear algebra, we know that Hd 2 — KerP c © ImP c . Therefore 
for every |a) G / Hd 2 , we have \a) = \a c ) + \a c >) where \a c ) G ImP c and \a c >) G KerP c with 
dim(ImP c )=(ii and dim(KerP c )=rf 2 — d\. Now we introduce the following Hermitian operator 

W = I dl ® I c - di|Vc)^c| - di(US ® I c M)(^ c p(LV ® I c ) (4.70) 

where I c = \i c ) (ic\ is the identity operator in the projected subspace of H^ 2 corresponding 

to a combination denoted by c. So the identity operator in the corresponding d\ <g> c?i subspace 
of c?i ® d 2 system is given as 

di-l 

l dl ®I c =J2\i,jc)(i,jc\ (4.71) 

i,j=0 

and l^c) is the the maximal entangled Bell-state in that subspace 

di-l 

l^) = ^=^K,i c ), (4.72) 

V "1 • n 
v 1=0 

in which |i c ) G ImP c for i = 0, • • • , d± — 1. It should be noted that the skew-symmetric 
operator is defined on the di-dimensional party of the d\ ® di system. To clarify, we give 
an example in 4 ® 5. Sine C\ — 5 then we have five 4-dimensional projected subspace for 
H 5 . The first one spanned by {|0i) = |0), |li) = |1), \2i) = |2), |3i) = |3)}, the second one by 
{|0 2 > = |0>, |1 2 ) = |1), |2 2 ) = |2>, |3 2 ) = |4)}, the third one by {|0 3 > = |0>, |1 3 ) = |2), |2 3 ) = 
|3), |3 3 ) = |4)}, the fourth one by {|0 4 > = |0), |1 4 ) = |1), |2 4 ) = |3), |3 4 ) = |4)} and the fifth 
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one by { | O5) = |1), 1 15} = |2), |2 5 ) = |3), |3s) = |4)}. Therefore one has five ways to construct 
the operator (14. 70 p . It is clear that the operator (14.701) is the same operator (12.31) in d\ <8> d\ 
bipartite quantum system which has been embedded in d\ ® d 2 Hilbert space. Consequently, 
the number of such operators obtained by embedding in this way is, in fact, ^s^' 

By the same arguments proposed in section 2, the operator (I4.70P can be considered as an 
entanglement witness which is defined on the d\ ® d 2 Hilbert space. Using instead of U d 
in (14. 70 p . we obtain a canonical form for W as 

W c = l dl ® I c - dil^cWcl - ® IcMc)(il>c\ TA (Jdi ® Ic), (4.73) 

where Qdi is an orthogonal matrix (QdiQj = Q J Qdi = I dl ) and J di is the canonical form 
of U di i.e., U di = Q di J di Qj. The action of (whose rank is 2n{) on the basis of the tri- 
dimensional single party Hilbert space is similar to (12. 16j) and (12. 17ft . Consequently, after some 
calculations, the expanded form of Wc becomes 

n\— 1 

W C =J2 d 2 *' ( 2 ^) c >< 2 *> ( 2 ^)J + l 2 *> ( 2 i + X ) c )( 2 ^ & + 1 U + I 2 * + 1> (2j)c)(2i + 1, (2j)J 

+ |2z + 1, (2j + l) c )(2z + 1, (2j + 1)J - \2i, (2z) c )(2j, (2j)J - \2i, (2i) e )(2j + 1, (2j + 1)J 
-|2i + 1, (2i + l) c >(2j, (2j) c | - \2i + 1, (2i + l) c )(2j + 1, (2j + 1)J 
-\2i + 1, (2j) c )(2j + 1, (2i)J + |2i + 1, (2j + l) c )(2j, (2z)J 
+ |2z, (2j) c >(2j + 1, (2i + 1)J - |2i, (2j + l) c )(2j, (2i + 1)J) 

711 — 1 dl — 1 

+ EE (l 2 ^ic)(2i,J c | + \2i + l,j c )(2i + 1,H + |j, (2i) c }(i ; (2*)J + \j, (2i + l) o )(j, (2i + 1)J 

-\2i, (2i) c )(j,j c \ - \2i + 1, (2i + l) c )(j,j c \ - \j,j c )(2i, (2i) c \ - \j,j c )(2i + 1, (2i + 1)J) 

di-l 

+ ^2(\hjc)(i,j c \-\hh)(j,j c \)- (4-74) 

These witnesses are optimal d-EWs for di < 4; otherwise, they are nd-EWs, so they can 
detect the entanglement of some PPT states. Optimality and nd-optimality problem for these 
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witnesses is similar to those ones which were discussed in subsection (2.3). Now we introduce 
PPT states in which some of them are entangled. Consider the following class of states 

j ni-l 

p = -(a \ip c )(ip c \ + aoY, (|2«, (2i) c )(2i, (2i) c \ + |2i + 1, (2i + l) c )(2i + 1, (2i + 1)J 

i=0 

-\2i, (2i) c )(2i + 1, (2^ + 1)J - \2i + 1, (2i + l) c )(2i, (2i) c \) 

n\ — 1 

+ ^ K+2,(2* )c \2i + 2, (2i) c )(2i + 2, (2i) c | + a 2m , (2l + 3 ) c |2i + 1, (2i + 3) c )(2z + 1, (2i + 3)J 

i=0 

-d(\2i + 2, (2i) c )(2i + 1, (22 + 3)J + \2i + 1, (2z + 3) c )(2i + 2, (2i)J)) 
m— l 

+ E «2i,( 2j )j2i, (2j) c )<2i,(2j) c | 

v£j=0,i—j^l,j-i^m-l 
ni— 1 

+ E a 2i+1 , (2j+1)c \2i + 1, (2j + l) c )(2i + 1, (2j + 1)J 

i^j=0,j-z^l,i-jyTii-l 

m— 1 

+ EK, (2j - +1) J2z, (2j + l) c )(2i, (2j + 1)J + a 2i+lj(2j)c |2i + 1, (2j) c )(2z + 1, (2j)J) 

i,j=0 

d\— 1 m— 1 

+ £ E«2^ c (|2j,0(2j;z;|+a 2j+M j2j + l,z;)(2j + l,i c | 

j=2ni j=0 

+ai,( 2j)c K, (2j) c )(i, (2j)J + a a2j+1)c |2, (2j + l) c )(i, (2j + 1)J) 

di-l 

i^jr'=2ni 
di— 1 d2 — di — 1 

+ E E a i,3 c >\h3c'){h3c'l ( 4 - 75 ) 

it is clear that the number of such states is C^ 2 . The positivity and PPT conditions are the 
same as for d <S> d cases except that we have an additional condition for positivity through the 
inequalities a^j t > (i — 0, • • • , di — 1 , j — 0, • • • , d 2 — di — 1). The subscript c' denotes 

c 

those states which lie in the KerP c (i — 0, • • • , d 2 — c?i — 1). Now the expectation value of the 
canonical EW with respect to the PPT state p gives out the following lower bound 

— 2r?i 

Tr(W cP )> 1 -, m = 2,3,4,.-.. (4.76) 

di + An\ + (d 1 — 2n 1 )(d 1 + 2ni — 1) 
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Therefore if those PPT states satisfy the following inequality 
— 2r?i 

dl + 4 „ f + (dl - 2 : i)(dl + 2 „ 1 - 1) ^^)<°. «. = 2.3,4,.., (4.77) 

then they are entangled. And, as before, if the states p violate the PPT conditions, then they 
always fulfill the following inequality 

4m (m - l) + (di - 2m) (di + 2m - 1) < % < ~ 2n i 



di + 2m ~~ di + 4n? + (di-2m)(di + 2m-l)' 

m = 2,3,4,-.. . (4.78) 

In the end of this paper, the idea of canonical EW introduced in equation (12 .lip can be 
extended to construct other EWs. Let us , at first, assume that the single party Hilbert spaces 
Hi and H 2 and the tensor product Hilbert space Hi © H 2 are defined on real field. As an 
illustration to this restriction, any entangled state which lie in a real tensor product Hilbert 
space Hi ®H 2 can be generated, by interactions or any entanglement generating process, from 
single party states which lie in real single party Hilbert spaces Hi and H 2 . It should be noted 
that any LOCC (local operation and classical communication) must be restricted on the real 
field. Therefore by these considerations, when we deal to detect the entanglement of a state 
which lies in a real tensor product Hilbert space H\®H 2 by constructing an EW, it is sufficient 
that our EW should have positive expectation value with respect to all real separable states. 
Now consider the following operator 

w = i d ® i d - d|vXVI - d(J T ® WXVPP ® id) 

- d(j' T (g) Lj^XVf^j' ® ^) - d(f T ® I d M)(i>\ TA (f ® l d ), (4.79) 



where 



J=j©j©j©... 
J'=j'©j'©j'©... (4-80) 
J =] ©j ©j ©..., 
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J 



> J 



(4.81) 



, in which the J, j' and j" are real skew-symmetric matrices in d- dimensional Hilbert space H 
with 

/ n i n n \ / q 1 \ / 1 

10 
0-100 
\ -1 / 

Clearly the ranks of J, j' and j" are An, n = 1, 2, 3, • • • with d = An + m where < m < 3. 
The action of J, j' and j" on the basis states of the d-dimensional single party subsystem are 
as 



10 
-10 
1 
\ 0-10/ 



-1 
-10 
\ 10 / 



J|4fc> = -\Ak + 1) j'\Ak) = -|4fc + 3) 

J|4fc + 1) = |4fc) j'|4fc + 1) = -\Ak + 2) 
J|4Aj + 2) = -|4fc + 3) j'|4fc + 2) = \Ak + 1) 



j"|4Jfe) = -|4fc + 2) 
j"|4Jfe + 1) = |4A; + 3) 
j"|4A; + 2) = |4/c) 



(4.82) 



J|4A; + 3) = |4A; + 2) f\Ak + 3) = \Ak) f\Ak + 3) = -\Ak + 1) 

The expectation values of the operator W with respect to the product states are as 



<r/l ® <CiW|77> <S) JO = 1 - | <Cl^7*> r - I <Cl J 



i<ciJ'k>r-i(ciJ"k>p 



(4.83) 



It is easy to see that the states \r)*}, J\rj), ^ \rf) and J"|?y) are orthogonal to each other when they 
are belong to a real single particle Hilbert space. Therefore, the expectation values become 
positive with respect to all real separable states. The relation between the operator (12. lip and 

<MB is 

W c = W + Dl A + Dl A , (4.84) 

where D^ 4 = d(j' T ®I d )|^)(^F A (J'®hi) and D^ 4 = d{f T ®l d )\^)(il)\ TA {f®l d ). The operators 
Dx and D 2 are positive operators so the PPT entanglement detection power of W is greater 
than Wc- Consider, for example, the following state with d = An which is defined in a real 
tensor product Hilbert space Hi ®%2 (and keeping in mind that any LOCC must be restricted 
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on the real field) as 

n-l 

= a c?|^)<^| + a ^[3(|4i,4i)(4i,42| + \Ai + 1, Ai + l)(Ai + 1, 4i + 1| 

i=0 

+ \4i + 2,4i + 2}(4i + 2,4i + 2\ + \Ai + 3, 4i + 3) (4i + 3, 4* + 3|) 
-|4i,4i)(4i + l,4i + l| - \Ai, 4i)(4i + 2, 4z + 2| - |4i, 4i)(4i + 3, 4z + 3| 
-|4i + 1, Ai + 1) (4i, Ai\ - \Ai + 1, 4i + l)(4i + 2, 4z + 2| - |4i + 1, 4i + l)(4i + 3, 4i + 3| 
-|4i + 2,4i + 2)(4i,4z| - \Ai + 2,4i + 2)(4i + 1, 4i + 1| - |4i + 2, 4i + 2)(4i + 3, Ai + 3| 
-|4i + 3,4z + 3)(4i,4i| - |4i + 3,4i + 3)(4i + l,4i + l| - \Ai + 3, 4z + 3)(4i + 2, 4i + 2|] 

n— 1 

+ ( a 4i,4j|4i,4j')(4i,4j| + a 4iAj+1 \Ai,Aj + l)(Ai,Aj + 1| + a 4iAj+2 \Ai,Aj + 2)(Ai,Aj + 2\ 

+a4i,4j+3|4i, 4j+3><4i, 4^+31+0^+1,4,-141+1, 4j)<4i+l, 4j|+a4i+i,4,-+i|4i+l, 4j+l><4z+l, 4j+l| 
+a 4 m,4j+2|4i + l,4j + 2)(4i + l,4j + 2| + a 4m , 4j+3 |4i + l,4j + 3)(4i + l,4j + 3| 

+a 4 *+ 2 ,4jl 4 * + 2, 4j)(4i + 2, 4 j | + a 4i+2Aj+ i\Ai + 2, 4j + l)(Ai + 2, 4j + 1| 
+04i+2,-y+2|4i + 2,4j + 2)(4i + 2,4j + 2| + a 4i+2 , 4j+3 |4z + 2,4j + 3)(4i + 2,4j + 3| 

+a 4 *+3,4,l 4 * + 3. 4 J')( 4i + 3. 4 J'I + a4i+3,4 i+ i|4i + 3, Aj + l)(4i + 3, Aj + 1| 
+a 4l+3 , 4i+2 |4z + 3,4j + 2)(4i + 3,4j + 2| + a 4l+3Aj+3 \4i + 3,4j + 3)(4i + 3,4j + 3|) 

rt— 1 

(«4i, 4 i+i 1 4i , 4i + 1 ) (4i , 4i + 1 1 + a 4i>4i+2 1 4i , 4i + 2) (4i , 4z + 2 1 
+a 4 i, 4 i+ 3 |4z,4i + 3)(4i,4i + 3| + a 4i+1Ai \Ai + 1, 4i)(4i + 1, Ai\ 
+a4i+i,4i+2|4i + l,4i + 2)(4i + l,4i + 2| + a 4i+Mi+3 |4i + l,4i + 3)(4i + l,4z + 3| 
+a4i+2,4i|4z + 2,4i)(4i + 2,4i| + a 4i+2 , 4m |4i + 2, 4i + 1) {Ai + 2, 4i + 1 1 
+a 4i+2)4i+3 |4z + 2, Ai + 3) (4i + 2, 4i + 3| + a 4i+3Ai \4i + 3,4i)(4i + 3,4i| 
+a 4 i+3, 4i +i |4i + 3,4i + l)<4i + 3,4i + 1| + a 4l+3Ai+2 1 At + 3 , Ai + 2) (4? + 3 , Ai + 2 1 ) 
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n-l 



+a ^(|4i + l,4z + 2)(4z + 3,4z| + |4i, 4z + 3)(4i + 2, 4z + 1| 

+ |4i + 3,4z)(4z + l,4i + 2| + |4i + 2, 4i + 1) (4z, 4z + 3|) (4.85) 
Its positivity and PPT conditions £1X6 £lS 

ao > 0, a4i+l,4t+2<i4i+3,4i > a 0' &4i,4*+3 a 4i+2,4i+l > ^0) 

z = 0,...,n-l (4.86) 

and 

0>M,4j a 4j,4i > a 0' a 4i,4j+l a 4j+l,4i > a (j' &4i,4j+2 a 4jr'+2,4i > Oq, 

a 4i,4jr'+3 a 4j+3,4i > ^q, ^4j+l,4j+l a 4j+l,4i+l > ^0) ^4i+l,4jr'+2^4j+2,4i+l > a 0! 

2 2 2 

a4i+l,4j+3«4j+3,4i+l > 0, Q , 04i+2,4j+2«4j+2,4i+2 > ^q, a4i + 2,4j+3&4j+3,4i+2 > Op, 

tt4i+3,4j'+3fl4j+3,4i+3 > a 0i 

i, j = 0, ...,n - 1, z^j, 

«4j+l,4ifl4i+3,4j+2 > «0' 04j,4i+lfl4i+2,4j+3 > a 0' * = 0' •"J 77 ' — ^! (4.87) 

respectively. The expectation value of W with respect to p, along with the positivity and PPT 
conditions, satisfies the following inequality 

Tr(Wp) > -4a n. (4.88) 

Clearly, by the relation ( 14.841) . the PPT entanglement detection power of Wc for this state is 
weaker than W. 



5 Conclusions 



In this paper, we have constructed new types of EWs, by using real skew-symmetric operators 
and maximal entangled state for bipartite d® d quantum systems analytically. The proving 
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of the positivity of expectation values of these witnesses has been done with respect to the 
separable states very easily. We have seen that by using various orthogonal transformations 
on canonical EWs one can obtain a large number of EWs. It has been shown that all of 
the canonical EWs in various ranks are optimal. When J is full-rank, the canonical EW is 
optimal nd-EW and for the cases which is not full-rank (2n < d) the corresponding canonical 
EW is not optimal nd-EW. In these cases, optimal nd-EW lies in the 2n ® 2n subspace of 
d® d Hilbert space. We have also constructed positive maps corresponding to canonical EWs 
by Jamiolkowski isomorphism. On the other hand, we have constructed the other types of 
witnesses in d ® d Hilbert space corresponding to the possible partitions of full-rank J. It has 
been shown that for a full-rank J there exist p(n) (the number of partitions of n) number of 
optimal EWs. Among these witnesses, there are one optimal d-EW and one optimal nd-EW. 
The other ones are composed of optimal nd-EWs and optimal d-EWs. We have also generalized 
our approach to the d\ ® d2 (d% < cfe) quantum systems. We have shown that there exist C d J 
distinct possibilities to construct EWs for a given d\ ® c?2 Hilbert space. In all of the cases 
where we have discussed, we have emphasized that the rank(J) > 4 (rank(Ji) > 4), because 
if we take rank(J) < 4 (rank(Ji) < 4) then we can not find out any nd-EWs. Also the idea of 
canonical EW has been extended to produce other EWs allowed to use only for detecting the 
entanglement of states which lie in a real tensor product Hilbert space. 

In each step, we have constructed a class of PPT states. The expectation values of 
constructed witnesses with respect to the corresponding PPT states give out a negative lower 
bound. This lower bound lies at the boundary of the entangled PPT states and entangled 
states that violate the PPT conditions. Finally we must mention that the other interesting 
issues remain unsolved such as for positive expectation values; we can not conclude that the 
corresponding states are separable or not. On the other hand, by replacing the W re d by the 
generalized reduction EW introduced in [27] , the generalized canonical EWs may be obtained. 
The work on these important points is under investigation by these authors. 
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Appendix A: 

Proving the lower bound ( 12.311) : 

The expectation value of the EW (12.221) with density matrix p (without normalization) in 
(I2.27P is given as 

n-l 

Tr(Wcp) = -(4n(n - 1) + (d - 2n)(d + 2n - l))a + ^ ( a 2i , 2j - + ^ 2i+h2j+1 + a 2i+lj2 . + a 2i>2j+1 ) 

n—l d—l n—1 d~l 

- 2 J2 Ci+ J2 + a ^ + a ^ + a ^+^> + o « ( A_1 ) 

i=0 i=2n j=0 i^j=2n 

Let us rewrite the the PPT condition in f l2.30p as 

a 2>, 2j a 2i,2i = $2i,2j a l 5 $2i,2j > 1 ; i, j = 0, n - 1 ; i ^ j 
a 2i+i,2j+i a 2j+i,2i+i = ^2i+i,2i+i a o i ^2i+i,2j+i > 1 ; i, J = 0, n — 1 ; i ^ j 
a 2l ,2 J+ i «2 J+ i, 2l = 5 2i,2j+i a o ; <Wj+i > 1 ; ij = 0, ...,n- 1 ; i^j 

a 2i+l,2i a 2i+2,2i+3 = ^2i,2i+2 a ' ^2i,2i+2«0 > Cj ; 1=0, .-.,n — 1 

a 2j,i a i,2j = $2j,i a o I 5 2 j,i > 1 ; j' = 0,...,n-l ; % = 2n, ...,d - 1 

a 2j+M fli,2i+i = £>j+i,i a o ! > 1 ; i = 0, n - 1 ; i = 2n, d - 1 

a»j %,i = ^ 2 j«o I s ij > 1 I i,j = 2n,...,d-l ; i^j (A-2) 

It is clear that 5 is invariant under permutation of its subscripts. It should be noted that 
for two variables with constant product, their summation becomes minimum when they are 
equal. The corresponding summation of each product in ( 1A-2I) appears in ( 1A-Ij) . Hence the 
minimum value for (1A-1|) is obtained when 



a 2i,2j = a 2 3 ,2i = $2i,2ja ; S 2i ,2j = <%,2i ; i,j = 0, n-l ; % ^ j 

a 2i+i,2j+i = a 2 3 +i,2i+i = ^2i+l,2j+lO0 ; <^2i+l,2j+l = ^2j+l,2i+l ! = 0, ...,U — 1 ; Z 7^ j 
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a 2«,2j+i = a 2j+i,2i — ^2i,2j+iao ; ^2i,2j+i = $y+i,2t ! «)J = 0, ...,n — 1 ; i ^ j 
a 2jjl = a »,2i = 8i,2jOo ; 5 ii2 j = 5 2 j,i ; J = 0, n - 1 ; i = 2n, d - 1 
a 2j+i,i = a «,2i+i = ^i,2i+i a o j "i,2j+i = <%+i,i ; J ' = 0, n — 1 ; % — 2n, d — 1 
Hi = a j,i = $i,j a o ; = ', i, j = 2n, d - 1 ; z ^ j 

Q = S 2ii 2i+2ao ; i = 0, n - 1 (A-3) 
Therefore the equations in (1A-3j) satisfy the PPT conditions and ensure to obtain minimum 



value for I I A- II) . Also the (1A-3I) along with the inequality Cj < 5 2 j+i, 21+3^0 = 0, n—1) gives 
the positivity conditions for p. To check this matter, one can write the first and second line of 
AA-3D for j = i + 1 and therefore a 2i+2 2i = ^,21+2^0 and a 2i+li2i+3 = <Wi,2i+3a Q (i = 0, n- 1) 
so a 2i+ 2,2i a 2i+i,2i+3 — ^2i,2i+2^2i+i,2i+30o > Cf (i = 0, ...,n — 1) which is the same as the positivity 
condition for p in ( 12.29ft . The equation ( 1A-1I) becomes 

n— 1 

Tr(W c p) = -(4n(n - 1) + (d - 2n)(d + 2n - l))a + a ^ (^, 2j + <W, 2j - +1 + <Wi, 2j + <W + J 

n—l d— 1 n—1 d— 1 

i=0 i=2n j=0 i^j=2n 

It is clear that we have used the invariancy of S with respect to permutation of its subscripts. 
Since Tr(Wc*p) is a linear strictly increasing function of various 5s then its minimum takes 
place in lower bounds of various 5s. Thus, when all of 5s are equal to one the lower bound of 
(IA-41) is (— 2nao). Consequently, by considering the normalization factor M (12.281) . which by 
(IA-31) gets its minimum value, we obtain 

— 2n 

^ W ^ d + tf + (d-2n)(d+2„-l) ; " = 2,3,4,.. (A-5) 

So we conclude that, the lower bound is obtained in the boundary of the entangled PPT states 
and entangled states that violate the PPT conditions. 

Appendix B: 
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Let us assume that H 2n <g> H 2n is a subspace of <g> Hilbert space. Wopc m (12.26ft . is a 
Hermitian operator in H 2n £S> H 2n which is embedded in the <g> and then its compact form 
is as 

Wopc = hn ® Ian - 2n\^){if>\ - 2n(J T ® I 2n ) |^> (^| Ta (J ® W (B-l) 

where 



Now we define |r/) <8> |C') G H2n <8> H 2n as the projection of products \rj) <g> |0 G <8> H d so the 



expectation value of Wopc with normalized products |ry) ® |C) is 

fal ® (CIWopcI^) ® |C) = (V'W)(C'\C) - KCV*)!* - \(CVW)\ 2 (B-3) 
By normalizing the states \r)') and |0)? we obtain the next equation 

(V\ O (CIWopcI^) ® |C) = (i/|i/>(CIO(i - KC'V*)I 2 - l(C"l%">l 2 ) (B-4) 

where 

w '^7m W) • lc " H vWy 10 (B " 5) 

The right hand side of the equation (1B-4I) . apart from the multiplier (VI 7 /) (010); i s similar to 
the right hand side of the equation (12.51) . Since J is full-rank in 2n- dimensional subspace then, 
by the same argument sketched in section 2, the above mentioned expectation value is zero by 
the products such as \n") ® (a\r)"*) + 0J\q")) and positive with other separable states so Wopc 
is an EW. To further illustrate, if \r)) ® \Q G H 2n <g> H 2n then the equation (1B-4|) becomes as 
the following equation 

(V\ ® moPclv) ® 10 = 1 - l(CI^)| 2 - l(C|J|^)| 2 (B-6) 

which was discussed in section 2. On the other hand, if \rj) <g> \Q G H d _ 2n eg) H d _ 2n , where 
Hd_2n ® H d _ 2n is the complement subspace of H 2 „ <g> H 2 „, then the expectation value of the 
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equation ( 1B-4I) with these products becomes zero because they have no component in H 2n ®H 2n . 
Consequently, the expectation value of Wopc with the following products is always zero, 

Pwopo = {|*7> ® 10 e H 2 „ ® H 2 „; |C) = a\ V *) + /?%>} 

U {\V) ® 10 e H d ® H d ; |C") = a|»H + /?,%")} U {|r/> <g> |C) G H,_ 2rt ® H d _ 2 „} (B-7) 
.where \rj") and |£") have been defined in equation ( 1B-5I) and |7/) <g> |£') is the projection of 



\r]) ® \Q in to subspace H 2n © H 2n and with other products is positive. 

Appendix C: 

As we considered in the paper, for a given partition fj, 2 , // 3 , • • • of n, the corre- 
sponding entanglement witness is given by 

W c (mi, A*2,M3, • • • ,/V> = W c (/ii,/i 2 ,/i 3 , • • ■ ,/i,)+0^(/i 1 )+0*- 4 (/i 1 ) + 0'- 4 (/i 3 ) + - ■ ■ +0*^-1) 
in which 

Wc(Mi, Ate, A* 3 , • • • , A^) = W c (/ii) © W c (/J 2 ) © W c (/i 3 ) © • • • © Wcr(/x„) 

Each of Wc(/i e ) (e = 1, 2, • • • , i>) is canonical EW on the 2/z e © 2/i e subspace of d ® d Hilbert 
space such as 



Wc(ft) = (|2z, 2j)<2«, 2j| + |2z, 2j + l)(2i, 2j + l\ + \2i + 1, 2j) (2i + l, 2j\ 

+ \2i + 1, 2j + 1) (2z + 1, 2j + 1| - |2z, 2z) (2j, 2j| - |2z, 2z) (2j + 1, 2j + 1| 
-|2i + 1, 2z + l)(2j, 2j| - |2z + 1, 2z + l)(2j + 1, 2j + 1| 
-|2i + 1, 2j){2j + 1, 2i\ + \2i + 1, 2j + l)(2j, 2i| 
+ |2z,2j)(2j + l,2i + 1| - |2i,2j + l)(2j, 2i + l\) 
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As shown in the paper, these witnesses for fj, e — 1 are zero and for the others are optimal 
nd-EWs. Consequently the Wc(fJ>i, 1^2, A*3> ' ' ' > /O is zero or optimal nd-EW and the operator 
0' A (/i e ) for (e = 1, 2, 3, • • • v — 1) is written as below 

hM E -l | _1 

0^(/i e ) = ^ ^ (|2i,2j)(2i,2j| + |2i,2j+l>(2i,2j + l| + |2i+l,2j)(2i+l,2j| 

i=Ml+M2H t-A J - e „ ! 3 = M1+M2H HMe 

+ |2z + l,2j + l)(22 + l,2j + l| -|2i,2i)(2j,2j| - |2i,2i)(2j + l,2j + 1| 
-\2% + 1, 2i + 1) (2j, 2j| - |2i + 1, 2i + l)(2j + 1, 2j + 1| 
|2j, 2i)(2j, 2*| + \2j, 2i + l)(2j, 2i + 1| + \2j + 1, 2i)(2j + 1, 2i| 
+\2j + l,2i + l)(2j + l,2i+l\-\2j,2j)(2i,2i\-\2j,2j)(2i + l,2i + l\ 
-\2j + l,2j + l)(2i, 2i\ - \2j + 1, 2j + l)(2i + 1, 2i + 1|) 
It is clear that the 0(/x e ) is positive operator so the 0* A ( / u e ) is optimal d-EW. 

Appendix D: 

In this appendix we give the positivity and PPT conditions for the following PPT state 
which corresponds to the partition fx 2 , // 3 , • • • , jjL v ) of n — |. 

where 

g(fii,H2, fis,'-' , AO = p(a»i) © p(^) © p(a»s) © • • • © pM 

Each of the p(fJ, e ) with (e = 1, 2, 3, • • • ,v) and cr(/u e ) with (e = 1, 2, 3, • • • , v — 1) is given as 
below 
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+a (|2i,2i)(2i,2i| + |2i + l,2i + l)(2i + l,2i + l| 

1=^+^2 +/u 3 H b" e _i 

-|2i, 2z)(2z + 1, 2i + 1| - |2i + 1, 2i + l)(2i, 2i\) 

/ij+Zijj+ZigH h/i e -l 

+ ]T (a2i+ 2) 2i|2i + 2,2i)(2i + 2,2i| +a2i+i )2 i+3|2i + l,2i + 3)(2i + l,2i + 3| 

i=/ij+/i 2 +/i 3 H h/i e _! 

— C<(j2i + 2, 2z)(2z + 1, 2i + 3| + |2i + 1, 2i + 3) (2i + 2, 2i|)) 

^j^+Zij+^gH hM E — 1 

+ J2 a 2i , 2j \2i,2j)(2i,2j\ 

i¥=3=Hi +M 2 +M 3 H l-M e -i -1 

^ +/i 2 +/i 3 H hM e -l 

+ £ a 2m>2j+ i|2i + l,2j + l)(2i + l,2j + l| 

*7^i=Mi +M 2 +M 3 H hM e _i ,J-»^lj»-jVA*e -1 

/ij+Zij+figH ^Me -1 

+ ^ (o2i,^-+i|2i, 2j + l)(2z, 2j + 1| + a 2i+1)2j |2z + 1, 2j)(2i + 1, 2j\) 

M=Ml+jU 2 +^ l-M e _i 

in which 



is the maximal entangled state in 2ji e <g> 2/j e subspace. and 

/i 1 +/i 2 +/i 3 +-+/i e -i f-i 

£ (a 2 . 2j \2i,2j)(2i,2j\ + a 2 . 2j+1 \2i,2j + l)(2i,2j + 1\ 

»=M 1 +M 2 +M 3 H ^M e _ 1 j = M 1 +M 2 +'^ 

+a 2i+1>2j |2i + 1, 2j)(2i + l,2j\ + a 2i+h2j+1 \2i + l, 2j + l)(2i + 1, 2j + 1| 
+a (|2i, 2^) (2j, 2j| + |2i, 2i) (2j + 1, 2j + 1| 
+ |2z + 1, 2% + 1) (2j, 2j| + \2i + l,2i+ 1) (2j + 1, 2j + 1|) 
+a 2 . 2i |2j, 2i)(2j, 2i| + a 2 . 2i+1 |2j, 2z + l)(2j, 2% + 1| 
+a 2 , +lj2i |2j + l,2z)(2j + 1,2*| + a m+1 \2j + l,2z + l)(2j + l,2i+ 1| 
+a (|2j, 2j)<2i, 2i| + \2j, 2j)(2i + 1, 2% + 1| 
+ |2j + 1, 2j + 1) (2i, 2*| + |2j + 1, 2j + 1) (22 + 1, 2i + 1|)) 
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where 



= ^A/"(aO + $>(aO 

e=l e=l 

with 

/i 1 +£t 2 +/j 3 H hMe" 1 A t i+A t 2+M 3 H l-Me^ 1 

N{n<) = 4a // e + K, 2 ,- + Ow+i^+i) + Yl ( a «,w+i + a 2i+ i, 

»^'=A»i+M2+M 3 H b« e _l *:i=^ 1 +M2+M3+"'+^e-l 

and 

/l^ +// 2 +// 3 H h/i e — 1 2 — 1 

»oo E E 

l°2i,2j + 2i,2j+l + a 2i+l,2j + a 2i+l,2j+l 

»=Mi +M 2 +M 3 +---+M e _i i=Mi +M 2 +M 3 H t"M £ 

~l" a 2j,2i + a 2j,2i+l + a 2j + l,2i + °2j + l,2i+l) 

The Positivity Conditions: 

The following conditions must be satisfied for each e — 1, 2, 3, • • • , v 
1) 

o > 

2) 

« 2i , 2j > f or i,j = Hi + H2^ \- fi e -i, - • - , Hi + fi2 -\ h/i e -l 

with i^j ; i - j ^ 1 ; j - i ^ /x e - 1 

3) 

a 2l+1 , 2j+1 > for i,j = + h A» e -i, • • • , A*i + A*2 H h /i e - 1 

wit/i i 7^ j ; j-i^l ; z - j ^ ^ e - 1 

4) 

a w l!y +i ^ ; a 2 i+ i,2j > f° r h j = A*i + A*2 H h /i e -i, • • • , /xi + //2 + • • • + He - 1 

5) 

a 2 l+2 ,2 l a 2 !+ i, 2l+ 3 > f or i = + H 1" • • • , A*i + A*2 H h (i e - 1 
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6) 

a ^ > ; a 2i,2 j+ i > ; a 2 i+ i,2,- > ; a 2i +i, 2 ,+i > 

a 2 ,, 2l > ; « 2 ,-, 2i+ i > ; « 2j+ i, 2l > ; ^ 

/or z = /ii+// 2 H hyu e _i, • • • ,^i+/i 2 H h/x £ -l and j = //i+/i 2 H h// e , • • • , ^-1 

Therefore these inequalities ensure that p(/xi, /x 2 , A*3> • • • , A*f) an d p(/-0 (e = 1, • are positive 
operators. 

The PPT Conditions: 

1) 

a 2i,2j a 2j,2i — °0 ! fl 2i+l,2j+l fl 2i+l,2i+l — tt ' a 2i,2j+l fl 2j+l,2> — O 

/or i,j = ^i+/i 2 H h // e -i,--- , Mi + M2 H h// e -l wit/i i^j 

2) 

a 2l+ i, 2l a 2 !+2 , 2l+ 3 > / or * = + H h • • • , Mi + ^2 H h /x e - 1 

3) 

2 . 2 . 2 , 2 

a 2i,2j a 2j,2i — % ! a 2 i, 2 j + i a 2j + l,2i — °0 ' a 2i+l,2j a 2j,2i+l — a ! a 2 i+l,2j + l a 2 j + l,2i+l — a 

/or 1 = ^+^ h^ e _i, • • • ,Mi+M2H h/i e -l and i = A*i+A*2H b«o " ' ' 

And these inequalities give the PPT conditions for /i 2 , jtt 3 , • • • , /x„) and p(/x e ) (e = 1, z/). 
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